Stoner-Wohlfart model applied to bicrystal magnetoresistance hysteresis 
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We calculate numerically the magnetization direction as function of magnetic field in the Stoner- 
Wohlfart theory and are able to reproduce the shape of the low-field magnetoresistance hysteresis 
observed in manganite grain boundary junctions. Moreover, we show that it is necessary to include 
biaxial magnetocrystalline anisotropy to fully describe the grain boundary magnetoresistance in 
La .7Sr .3MnO3/SrTiO3 bicrystal tunnel junctions. 
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Can the fairly simple model presented by Stoner and 
Wohlfart in 1948 1] provide insights in the processes that 
occur in manganite bicrystal grain boundary junctions? 
In the following we show how the coherent rotation of 
the magnetization vector can explain the characteristic 
features in the field dependence of the magnetoresistance. 

Recently, a substantial low- field magnetoresistance was 
observed in grain boundaries of perovskite manganites 
0,13 and other half- metallic ferromagnets 0,0- In order 
to understand and exploit this effect several studies have 
used bicrystal grain boundary junctions Since 
bicrystals are samples with well defined crystal orienta- 
tions with respect to an interface, they are in a way ideal 
systems for studies of the behaviour of magnetic tunnel 
junctions. Traditionally, the conductivity of magnetic 
tunnel junctions is compared with Julliere's model [9j, in 
which the tunneling magnetoresistance is determined by 
the spin polarization of the electrodes. A more realistic 
scenario was considered by Slonczewski [lfjj, who used a 
method to match the wave functions across the tunnel- 
ing barrier. He derived an expression for direct tunneling 
through the interface which includes the angle between 
the directions of the magnetization of the electrodes. 

However, direct tunneling is not sufficient to explain 
the measured transport data in manganite grain bound- 
ary junctions. In fact, even though several scenarios have 
been suggested to explain measured data, the transport 
mechanism in bicrystal grain boundaries of manganites is 
still not fully understood. Klein, Hofener and coworkers 
[ill pointed out the strong impact of inelastic pro- 
cesses in the barrier region. They concluded that multi- 
step tunneling via a number of localized states within the 
barrier had to be added to the elastic tunneling conduc- 
tivity contribution. Several different inelastic processes 
have been suggested in order to explain experimental 
data, most of them emphasizing scattering at magnetic 
intra-barrier states Il3l llj, lla . 



So far, most studies of the transport mechanism in 
manganite grain boundaries have focused on the non- 
linearity of the / — 1^-curves, the magnitude of the low- 
field magnetoresistance or the shape of the magnetoresis- 
tance at high fields. The shape of the low-field hysteresis, 
directly related to the magnetization reversal process at 
the grain boundary, has long been neglected. Recently 
however, one attempt to reproduce the hysteresis was 
presented: Garcia and Alascio [16j minimized the to- 
tal magnetic energy and obtained the magnetoresistance 
from the magnetization directions. They considered the 
energy within the grain boundary region, and assumed 
uniaxial anisotropy in order to reproduce the experimen- 
tal data. 

In this paper we explore the possibility to describe 
the low-field magnetoresistance in a bicrystal junction 
in terms of coherent rotation of the magnetization direc- 
tions, similar to what was presented more than half a 
century ago by Stoner and Wohlfart In a numerical 
calculation we demonstrate that the shape of the hys- 
teresis in the magnetoresistance can be obtained by a 
simple energy-minimization technique, where each grain 
is treated individually. We show that a biaxial magne- 
tocrystalline anisotropy very well reproduces the mea- 
sured magnetoresistance curve in a bicrystal grain bound- 
ary junction. 

In the original paper by Stoner and Wohlfart from 
1948 Q, they described the magnetization reversal of 
uniformly magnetized ellipsoidal particles. Generally the 
equilibrium domain structure and magnetization reversal 
processes are determined by the balance between the ex- 
change, magnetocrystalline, magnetostatic and Zeeman 
energy contributions (see e.g. Ref. |l7f). Let us con- 
sider the two-dimensional case of a microbridge crossing 
a single bicrystal grain boundary, a type of device that 
has been well studied experimentally. The easy axis of 
magnetization, the magnetic field B and the magnetiza- 
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FIG. 1: The directions of the easy axis, the applied magnetic 
field and the magnetization (a, j3 and 7, respectively) are all 
defined with respect to the grain boundary. The magnetore- 
sistance is measured with the current along the length of the 
structure. 



tion M are oriented in directions (a, (3 and 7, respec- 
tively) defined relative to the grain boundary as depicted 
in Fig.^ Furthermore, let us first consider the two single 
crystal electrode sides (left and right) as being uniformly 
magnetized and decoupled from each other, i.e. only the 
variation of the magnetocrystalline and Zccman energies 
have to be considered. Hence the field dependent part of 
the energy density for one of the sides can be written as 



K 

uj = — sin 2 2(a - 7) - MBcos(/3 - 7), 



(1) 



where K is the (first order) biaxial anisotropy coefficient. 
In the case of uniaxial anisotropy the first term should 
be replaced by 0] 

K u sin 2 (a - 7). 

For fixed values of K and M and the angles a and 
P, the magnetization direction 7 can be obtained from 
Eq. Hby tracing a local energy minimum in the uj(B,"f) 
diagram. The starting point in this numerical method is 
chosen to be a global minimum at high magnetic field. By 
stepping the magnetic field we obtain "f(B), as shown in 
Fig. |2 With the assumption that magnetization reversal 
occurs by coherent rotation and \M\ = M s , where M s is 
the saturation magnetization, we can trace the rotation 
of M and find M(B). 

To apply this model to the bicrystal grain boundaries 
we trace the local minima defined by Eq. ^ f° r the left 
and right side independently. Hence we obtain jl(B) and 
Jr(B) for the two sides with easy axis in the direction 
oth and an, respectively. In Fig. [3|we show 7l(-B) and 
7fl(B) calculated for biaxial anisotropy. 

Now, for comparison with transport data, we con- 
sider the two sides connected by a tunneling barrier 
at the grain boundary. This yields a ferromagnet/non- 
magnet/ferromagnet (F/N/F) spin valve structure. As 
theoretically demonstrated by Slonczewski , the spin- 
dependent tunneling at the grain boundary is propor- 
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FIG. 2: Constant energy contours (u>) calculated from Eq. 
for constant K and M (KjM = 38 mT) using a) uniaxial 
and b) biaxial anisotropy, with a = 60° and j3 = 90°. The 
bold line indicates a local minimum traced from —B max to 
+B m «, and the direction of rotation is indicated by arrows. 



tional to the cosine of the angle between the magnetiza- 
tion directions of the electrodes. Thus the spin dependent 
conductivity is 

(J sp = G sp [l + P 2 cos(7 L - 7^)], 

where G sp is the spin polarized conductivity at jL~"fB, — 
it/2, and P is the spin polarization of the electrodes. 
With a non-spinpolarized contribution a ns = G ns , the 
total resistivity of the tunnel barrier is 
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P : 



1 + P 2 cos(7 L - 7fl ) + G' 



(2) 



where G = G ns /G sp . Spin- flipping inelastic processes 



can also be included in <j r . 



For each value of B we 
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FIG. 3: The angle of magnetization 7 for the left and right 
electrodes with or — —oil = 60°, (3 = 90°, and K/M = 
38 mT calculated for the case with biaxial anisotropy. The 
magnetic field starts from a high negative field and sweeps to 
positive values. 
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FIG. 4: The magnetoresistance of a manganite grain bound- 
ary. The circles are measured at 1.9 K according to the de- 
scriptions in the text. The solid line is a fit to Eq. with 
P = 1 and G = 20 (using K/M = 38 mT in Eq. [TJ. The 
arrows indicate the shoulder discussed in the text. 



calculate jl and 7^ and hence obtain the resistance hys- 
teresis p(B). 

The results of the numerical simulation are compared 
to the measured magnetoresistance of a bicrystal sample. 
The sample is a Lao.7Sro.3Mn03 (LSMO) film grown on a 
SrTiC>3 bicrystal substrate. The bicrystal is a symmetric 
(OOl)-tilt grain boundary with a misorientation angle of 
±15°, i.e. the [100] directions of the left and right sides 
are in-plane rotated 15° (in opposite directions) with re- 
spect to the grain boundary. The 120 nm thick film was 
grown by pulsed injection metal-organic chemical vapour 
deposition (MOCVD). More details on the growth proce- 
dure can be found in Ref. The good epitaxy of the 
film was verified by x-ray diffraction in 9-28- and c/>-scans. 
The out-of-plane lattice parameter of the film was found 
to be 3.86 A, to be compared to the bulk pseudo-cubic 
lattice parameter of the perovskite structured LSMO of 
3.88 A. It is known that this kind of tensile strain leads to 
an in-plane biaxial magnetocrystalline anisotropy in the 
(110)-directions.|Il Hence, this LSMO on SrTi0 3 
system can be analyzed with the two-dimensional model 
considered for the Stoner-Wohlfart theory. The satura- 
tion magnetization of the film, determined in a super- 
conducting quantum interference device (SQUID) mag- 
netometer, was 400 kAm -1 . 

A 5 fim wide and about 150 fim long microbridge cross- 
ing the grain boundary was defined by photolithography 
and Ar-ion milling, and the ends of the microbridge were 
connected to electrical contacts. The resistance of the 
grain boundary was measured in a four-contact geome- 
try with a current bias of 10 /iA. The magnetoresistance 
(R(B)/R(0)) was measured with the field applied along 
the length of the microbridge, i.e. in the plane of the 



film and perpendicular to the grain boundary. The mag- 
netic field was swept with 0.05 T/min. The dependence 
of the magnetoresistance on the magnetic field direction 
is presented in [22j- 

The two energy terms in Eq. ar e comparable in size 
when the field is B w K/M. For LSMO films on SrTiOs 
numerical values of K of 1.6 - 5.7xl0 3 J/m 3 [U Ei| 
have previously been reported. In addition we have M = 
400 kAm -1 , as stated above. Hence a simple estimate of 
the critical field for magnetization switching would give 
K/M =4 - 14 mT. 

Since the measured sample has a misorientation angle 
of ±15°, and the easy axes are in the [110] and equiva- 
lent in-plane directions, (Xl.r — ±60°. The field applied 
along the microbridge corresponds to = 90°. Using 
these parameters in Eq.fjwe can fairly well reproduce the 
measured grain boundary magnetoresistance, see Fig. 0] 
Here we use a K/M ratio of 38 mT, P = 1 and G = 20. 
Hence, the K/M from the experimental results is within 
one order of magnitude of the value estimated above. 

The width of the simulated magnetoresistance hystere- 
sis (the distance between the peaks) is uniquely deter- 
mined by the K/M ratio. The general shape of the 
curve is determined by the chosen type of magnetocrys- 
talline anisotropy (uniaxial or biaxial) together with the 
angles of the easy axes (a^ and an). The curvature of 
the ma gne toresistance mainly depends on the transport 
model |23j. The height of the simulated magnetoresis- 
tance depends on the values chosen for P and G. It is 
obvious from Eq. that a decreased spin polarization 
(P < 1) has virtually the same effect as an increased 
non-spinpolarized current (G > 0). In the simulation we 
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used P = 1 and G — 20, but for instance G = and 
P = 0.22 would yield the same magnetoresistance mag- 
nitude. Without further information about either P or 
G, a ratio between the two can not be extracted from the 
measured data with our model. 

We note that in the measured curve, there is a peak 
and a shoulder on each wing of the hysteretic curve. Both 
these features are represented in the curve (solid line in 
Fig. 0} calculated from Eq. [3 Each feature is associated 
with a jump in the magnetization direction. A closer look 
into the model of biaxial anisotropy reveals that the first 
jump comes from the switching from a high energy state 
to one with lower energy, marked I in Fig. In the 
low-energy valley there are two local minima. Thus the 
first part of the magnetization rotation is paused when 
the first of the two minima is reached in the direction 
of rotation. However, under some geometrical conditions 
there exists a state even lower in energy, and the second 
jump (II) is a transition to that state. The dwelling in 
the intermediate state creates the shoulder (marked by 
arrows in Fig. 0J) in the magnetoresistance curve. We 
note that this intermediate state does not exist in the 
case of uniaxial anisotropy (see Fig. [5^,). Thus, we con- 
clude that with the Stoner-Wohlfart model we can well 
reproduce the characteristic features of the low-field mag- 
netoresistance related to bicrystal grain boundaries, in 
the case when the magnetic field is perpendicular to the 
grain boundaries and a biaxial anisotropy is included. 

We find the agreement between the numerical values 
obtained in the model and in experiments to be good, al- 
though there are details that are not fully reproduced. In 
the model, the magnetocrystalline anisotropy K and the 
magnetization M uniquely determine the position of the 
peak. Considering the uncertainty in the experimental 
material parameters the estimated position agrees well 
with the experiments. By including surface effects at 
the grain boundary, which may locally decrease M and 
enhance the effective value of K, the agreement can be 
improved. 

Our simulation underestimates the curvature of the 
magnetoresistance curve. The curvature is determined 
by the transport model and here we have used the one 
by Slonczewski [IH ]. who considered a single-band direct 
tunneling process. From I — V measurements we know 
that direct tunneling is insufficient to explain the trans- 
port properties in this kind of junctions [Tl|. We also note 
that an enhanced curvature has been observed with an 
increased bias current [2||. Hence, our results propose 
that a more complete transport model for this kind of 
magnetic tunnel junctions should be developed. Further- 
more, the good agreement between the Stoner-Wohlfart 
model and the shape of the magnetoresistance curves mo- 
tivates future studies of the electrical transport in this 
kind of systems. 

The presence of the shoulder in the magnetoresistance 
hysteresis and its relation to the biaxial anisotropy have 



not, to our knowledge, been explicitly studied previously. 
Garcia and Alascio 16] did use a Stoner-Wohlfart-like ap- 
proach to the problem, but failed to demonstrate the in- 
fluence of biaxial anisotropy. However, Philipp et al [24| 
and Todd et al |25j present magnetoresistance curves for 
single magnetic bicrystal junctions with the field applied 
perpendicular to the grain boundary, and both curves 
show a behaviour similar to what can be expected from 
a coherent rotation of magnetisation direction with biax- 
ial magnetocrystalline anisotropy. 

A reason for the shortage of experimental data may be 
that the shoulder only appears under two conditions: at 
low temperatures and for a certain range of angular re- 
lations between the grain boundary, the anisotropy axes, 
and the ap plie d magnetic field. In the paper by Gunnars- 
son et al. |22j, which shows measurements on the same 
sample as presented here, the shoulder is absent in the 
magnetoresistance hysteresis measured at 100 K (with 
the field perpendicular to the grain boundary) . One rea- 
son for the T-dependence is of course the strong vari- 
ation of the magnetic properties with temperature. It 
should be noted that the Stoner-Wohlfart model is valid 
for coherent rotation and thus sufficient to describe the 
case when the field is applied perpendicular to the grain 
boundary, in our experimental setup. With the field 
applied in another direction (e.g. parallel to the grain 
boundary j2^|) the magnetization reversal processes be- 
come more complex. In that case one has to include con- 
tributions from the magnetostatic and exchange energy 
terms in the simulation. 

In summary, we proposed that the Stoner-Wohlfart 
theory can be applied to manganite grain boundary mag- 
netoresistance. Our conclusion is that the model of two 
magnetically decoupled electrodes with coherently rotat- 
ing magnetization vectors very well explain the experi- 
mental data. We find that in the case of Lao.7Sro.3Mn03 
junctions on SrTi03 bicrystals it is necessary to include a 
biaxial magnetocrystalline anisotropy to fully reproduce 
the measured data. In addition our study shows that 
a bicrystal grain boundary, with its well defined angles 
of magnetization, constitutes a system well designed to 
explore the physics of magnetic tunnel junctions. 
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